Fractional topological phases and broken time reversal symmetry in strained graphene 



Pouyan Ghaemi,^' ^' ^'j^] Jerome Cayssol,^'"^' ^ Donna N. Sheng,^ and Ashvin Vishwanath^' ^ 

^Department of Physics, University of Illinois, Urbana, IL 61801 
^Department of Physics, University of California, Berkeley, California 94720, USA 
^Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720 
^Univ. Bordeaux and CNRS, LOMA, UMR 5798, F-33400 Talence, France 
^ Max- Planck- Institut fiir Physik Komplexer Systeme, Nothnitzer Str. 38, 01187 Dresden, Cermany 
^Department of Physics and Astronomy, California State University, Northridge, California 91330, USA 

We show that strained or deformed honeycomb lattices are promising platforms to realize frac- 
tional topological quantum states in the absence of any magnetic field. The strained induced pseudo 
magnetic fields are oppositely oriented in the two valleys [iHS] and can be as large as 60-300 Tesla as 
reported in recent experiments [4] [5]. For strained graphene at neutrality, a spin or a valley polarized 
state is predicted depending on the value of the onsite Coulomb interaction. At fractional filling, the 
unscreened Coulomb interaction leads to a valley polarized Fractional Quantum Hall liquid which 
spontaneously breaks time reversal symmetry. Motivated by artificial graphene systems [5^1, we 
consider tuning the short range part of interactions, and demonstrate that exotic valley symmetric 
states, including a valley Fractional Topological Insulator and a spin triplet superconductor, can be 
stabilized by such interaction engineering. 



PACS numbers: 73.43.-f,73.22.Pr,72.80.Vp 

Fractional Quantum Hall (FQH) phases are macro- 
scopic scale manifestations of quantum phenomena with 
unique features including the fractional charge and 
statistics (abelian or nonabelian) of elementary excita- 
tions. This topological order originates from the strong 
Coulomb interactions between electrons moving in a par- 
tially filled Landau level induced by a strong magnetic 
field. Recently Chern insulator models with a nontriv- 
ial flat band [9HTT] were also shown to exhibit topolog- 
ical order in the absence of any magnetic field p!2UT7] . 
Those so-called Fractional Chern Insulators (FCIs) ex- 
plicitely break time-reversal symmetry T as did the orig- 
inal Haldane model [18]. In contrast, fractional topo- 
logical insulators (FTIs) p!9H22] can be naively thought 
of as two copies of time-reversed Laughlin FQH states, 
thereby obeying time reversal symmetry T. In spite of 
few proposals [23^, , the experimental implementation 
of FCIs and FTIs remains very challenging. 

Motivated by recent experimental advances [H |5] , we 
introduce another route towards fractional topological 
phases making use of the gauge fields that can be gen- 
erated in a deformed honeycomb lattice [1-3 . The as- 
sociated effective magnetic fields are opposite in the two 
different valleys and therefore they do not break the time 
reversal symmetry T [1 . Indeed, a scanning tunneling 
spectroscopy study [4 confirmed that straining graphene 
could yield flat Pseudo Landau Levels (PLLs) [2, 3 with 
effective fields as high as 300 T in each valley. Most 
recently by designing a molecular honeycomb grid of 
carbon monoxide molecules on top of a copper surface. 
Gomes at al. [5 were able to observe the linear dispersion 
of Dirac fermions in graphene, and furthermore to gener- 
ate nearly uniform pseudo-magnetic fields as high as 60 
T by deforming this grid [5] . Finally other realizations of 
artifical graphene systems, in patterned GaAs quantum 



wells [6] or with cold atoms trapped in hexagonal opti- 
cal lattices [7", "8 , also provide experimental platforms to 
create strong valley-dependent effective magnetic fields. 

In this Letter, we first consider real graphene under 
strong pseudo-magnetic fields generated by a mechanical 
strain. We investigate the interaction-driven phases in 
the n = PLL using mean field and numerical exact- 
diagonalization. The unscreened Coulomb interaction 
stabilizes a valley polarized Laughlin liquid at filling 2/3 
of the n = PLL. This states breaks spontaneously 
time reversal symmetry and is characterized by a finite 
charge Hall effect. At the neutrality point, we predict 
that strained graphene have either spin polarized or a 
valley polarized state, depending on the strength of the 
on-site interactions, with current estimates [25 favoring 
the former state. Second we have investigated what type 
of interactions could destabilize the valley-polarized 2/3 
state towards a valley-symmetric (time-reversal invari- 
ant) FTI. It turns out that the 2/3 state is rather robust 
for realistic interactions. Neverthess attractive local cor- 
rections to the Coulomb interaction can stabilize this val- 
ley FTI, which is a FTI where the valley plays the role 
of spin. Finally further increase of the attractive inter- 
actions leads to a spin triplet superconductor. Since the 
reported effective magnetic field strengths are around 300 
T [4 or 60 T [5 , the predicted phases might conceivably 
be realized with larger energy gaps than in FQH states 
under a real magnetic field. 

ModeL The noninteracting part of our model has 
been proposed by Guinea et al. in order to realize PPLs 
in strain graphene under zero magnetic field [2]. The 
corresponding tight-binding Hamiltonian reads 

H0 = Y. {t + 5ta{Vi)){0){Yi)b{Vi+5a)+h.C.), (1) 

a=l,2,3 
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where 6ta{Ti) is the strain-induced variation of the near- 
est neighbor hopping amphtude (with respect to the un- 
perturbed value t 2.7 eV) between A-site at and 
B-site at + 6a of the bipartite honeycomb lattice [1 . 
The smooth deformation field Sta{Ti) is chosen in such 
a way to produce a nearly uniform magnetic field with 
a valley-dependent sign [21 |26]. The valley dependent 
vector potential A^(r) = ^ X]ci=i 2 3 ^^a(r)e*^"^^ mini- 
mally couples to linearly dispersing low energy excita- 
tions near the Dirac points located at momenta with 
K = (47r/3>/3ao)ex and ^ = ±1, ao being the carbon- 
carbon bond length [1 . The uniform pseudomagnetic 
field induces a pseudo Landau level (PLL) electronic 
structure En = ^^y2ehvpB\n\, where n is the relative 
integer labelling the nearly flat levels (see supplemen- 
tary). Beside the macroscopic orbital degeneracy, each 
of those PLLs has a four-fold degeneracy associated with 
the spin and valley isospin degrees of freedom. In con- 
trast to the full SU{4:) symmetry of graphene in an ex- 
ternal real magnetic field [27l |28] , the internal symmetry 
of strained graphene is SU{2) for the spin and only Z2 
for the valley degree of freedom. 

In this work, we study interaction effects within the 
partially filled zero-energy flat band (n = PLL) cre- 
ated by strain. We consider the following interaction 
Hamiltonian on the honeycomb lattice: 

Hint= ^ V{Yi-Yj)n{Yi)n(Yj) ^UQ^n(Yi)n(Yi) 

(2) 

^Unnn ^ n{Yi)n(Yj), 

where F(r^ — r^) = e^/47re|ri — rj| denotes the bare 
Coulomb potential, n(r^) the fermion number operator 
on site r^, and (r^,rj) represents summation over all 
pairs of next-nearest-neighboring (NNN) sites. The bare 
Coulomb interaction is the dominant interaction due to 
the poor screening in neutral graphene. Nevertheless we 
also allow arbitrary modification of the short-distance 
part of the Coulomb interaction by adding local on-site 
and NNN interactions with respective strengths Uq and 
Unnn- The nearest neighbor interaction is not effective 
in presence of strong pseudo-magnetic field because in 
the zero energy PLL the noninteracting wave functions 
are localized on a single sublattice (see Supplementary). 
Interesting proposals for altering short ranged interac- 
tions using substrates with momentum dependent di- 
electric susceptibility has been discussed [29 . Unfortu- 
nately the actual values of /7o and Unnn are not known 
in strained graphene although first principles calculations 
yield total on-site coupling Uq = 9.3 eV and a small devi- 
ation Unnn — — 0.04e^/47reao of the NNN coupling from 
its bare Coulomb value in freestanding (and unstrained) 
graphene in zero magnetic field [25^. 
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FIG. 1: (Color online). The n = PLL at fractional 
filling factor 1/ = —2 + 2/3: spinless electrons. Up- 
per panel left: Low energy spectrum as a function of the 
next-nearest-neighboring coupling Unnn between opposite val- 
leys (deviation from the pure Coulomb value). In the region 
—0.73 < Unnn < —0.58, the nine lowest energy states be- 
come close together and almost degenerated, thereby forming 
the groundstate manifold (GSM) of the valley FTI. The inset 
shows the groundstate energies of the valley polarized (V- 
2/3 FQH) and FTI states. Upper panel right: The boundary 
phase dependence and the robust gap between the GSM and 
higher energy states for Unnn — —0.6 inside the FTI phase. 
Lower panel: Phase diagram as a function of Unnn for spin- 
less electrons. Parameters for the exact diagonalization: The 
noninteracting orbitals are determined on a 24 x 24 lattice 
with a pseudomagnetic flux $o/48 per hexagon (see supple- 
mentary material). The degeneracy of n = PLL is A^s = 12 
per spin direction and per valley. The low energy spectrum 
is calculated for Ne = S {Nl = Nr — 4) electrons with polar- 
ized spin occupying those Ns = 12 states. Energies are given 
in units of e^/47reao ^ 10 eV. 



Practionalized phases and superconductivity at 
2/3 filling of the n = PLL. Fractional Hall states in 
graphene under an external magnetic field were reported 
experimentally [3QH32] . Although strain produces flat 
PLLs, it is not evident that interactions can generate in- 
compressible phases at fractional filling in time-reversal 
invariant strained graphene. We focus here on the 2/3 
filling of the four- fold degenerate n = PLL. This 2/3 fill- 
ing has been studied so far in graphene sheets [33l [34] and 
in GaAs Hall bilayers [35] under real magnetic field. In 
the present case of strained graphene, this particular fill- 
ing allows for interesting possibilities including valley fer- 
romagnetism (which breaks spontaneously time-reversal 
symmetry), valley symmetric topological states, and also 
superconductivity. 
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FIG. 2: (Color online) The n = PLL at fractional filling 
factor v — —2 + 2/3: spinfull electrons. Left panel: The 
energy of different ground states as a function of the next- 
nearest-neighboring coupling Unnn defined in Eq. ([2]). In the 
range Unnn < —0.8, which includes the pure Coulomb inter- 
action of realistic graphene (Unnn = 0), the ground state is 
a valley polarized and spin singlet FQH state (green crosses) . 
Only a very significant attraction Unnn < —0.8 can destabi- 
lize this state towards a valley unpolarized and spin polar- 
ized superconducting state (red crosses). This superconduct- 
ing phase is characterized by a finite superfluid density as 
shown in the inset. Right panel: Two lowest energies in each 
momentum sector as a function of k for the valley-polarized 
state (Unnn = —0.6, upper) and for the spin-polarized su- 
perconductor {Unnn = —1, lowcr) . Parameters foT the exact 
diagonalization: same than in Fig. 1 but with the spin degree 
of freedom. 

Real graphene: time reversal breaking FQH state in a 
single valley. We first consider real graphene with the 
unscreened Coulomb interaction, namely Uq = Unnn = 
in Eq. ([2|. Then the ground state is found to be a val- 
ley polarized FQH state both for spinless (Fig. [T]) and 
spinfull (Fig. [2| electrons. This valley-polarized state 
breaks spontaneously the time-reversal invariance of the 
strained graphene Hamiltonian, and spins are in a singlet 
state as in the 2/3 FQH states [33H35] obtained under real 
magnetic field. Due to the large values of strain-induced 
pseudomagnetic fields, this state may be realized with 
elevated energy scales, allowing for the stabilization of 
fragile states. In order to test quantitatively the robust- 
ness of the 2/3 valley polarized FQH state, we now vary 
the parameter Unnn in the Hamiltonian Eq. ([2| . It turns 
out that the 2/3 valley polarized state is rather stable 
both in the spinless (Fig. [T]) and spinfull (Fig. [2| cases. 
Nevertheless when Unnn is sufficiently negative, exotic 
valley symmetric phases can also be realized as detailled 
below. For clarity we describe separately the spinless and 
spinfull cases. 

Spinless fermions and valley fractional topological in- 
sulator (FTI): Let us consider spinless electrons and de- 
compose the NNN coupling of Eq.([2| into an interaction 
between opposite- valley electrons {U^^^) and an interac- 
tion between same- valley electrons {U^^^). We first tune 
the intervalley correlations (U^nn) while U^^^ = (but 



note that electrons in the same valley still interact via 
the bare Coulomb potential). 

In some intermediate parameter range (—0.73 < 
Unnn < —0.58), an interesting quantum phase emerges 
with nine nearly degenerated states forming a ground 
state manifold GSM (Fig. [l]a, lines with symbols) which 
is well separated from the higher energy states (Figjlja, 
lines without symbols). This valley-symmetric and 9-fold 
degenerated phase is called here valley fractional topolog- 
ical insulator, since valley plays the role taken by spin in 
the previously discussed "spin" FTIs p!9H22] . Moreover 
the momentum quantum numbers of these states are at 
k = and other k determined by shifting the momen- 
tum of each electron by an integer multiple of 27r/A^s, 
where Ng is the PLL orbital degeneracy. This determines 
three different momenta sectors {k = 0, 7r/3, and 27r/3) 
and there are three near degenerate states in each sector. 
These sectors can be considered as ground state flows 
from one sector to another upon inserting flux through 
adding the twist boundary phase (Fig. [l]b). 

As a complementary characterization of the valley FTI 
phase, we further perform valley-pseudospin Chern num- 
ber calculation [36l [33 by adding the same boundary 
phase along x-direction, and the opposite ones along 
^/-direction for both valleys [SB] |39]. This generalized 
pseudospin Chern number is well defined as the electron 
number in each valley is conserved thus that the valley- 
pseudospin is a good quantum number. We find a total 
Chern number quantized to 6 for all nine levels, charac- 
terizing the 2/3 fractionalized valley spin-Hall effect. 

Finally we can also turn on and increase the intraval- 
ley part of the NNN coupling U^^^ (see supplemen- 
tary). In the limit of very large intravalley correla- 
tions {U^^^ oo), we expect a totally valley-decoupled 
1/3 + 1/3 phase consisting of two 1/3 Laughlin FQH 
states with opposite chiralities. We find no phase tran- 
sition between the valley FTI state discussed above (at 
^nnn — ^ ) dccouplcd fractional valley Hall in- 

sulator (see supplementary). 

The above results for spinless electrons can be summa- 
rized in a phase diagram (Fig. [l]c). For U^^^ > —0.35 
(which includes realistic graphene), electrons have a nat- 
ural tendency towards valley ferromagnetism, which is 
expected for repulsive interactions in a such a flat band 
system. In order to realize the valley FTI, one needs 
to counteract this trend by tuning an attraction between 
electrons in the opposite valleys. Besides, one also notices 
a narrow range of parameters (—0.58 < U^f^^ < —0.35) 
where the valley-polarization is lost but the GSM de- 
generacy not yet achieved. The understanding of this 
crossover region between the valley polarized FQH insu- 
lator and the valley FTI is still lacking and will be studied 
elsewhere. Finally superconductivity might appear when 
attraction is dominant {U^^^ < —0.58). This flat band 
superconductivity is discussed below in more details for 
the spinfull electrons. 
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Spinfull fermions and spin triplet superconducting 
state. We now consider spinfull fermions and we tune 
Unnn without distinguishing the valleys. For sufficiently 
large attraction (Fig. |2|, namely Unnn ^ —0.8 (note that 
when added to the Coulomb repulsion, this ends up giv- 
ing a somewhat smaller but still attractive next nearest 
neighbor interaction of /7*°*^ = —0.2), the ground state of 
the spinful model becomes a spin triplet and valley singlet 
superconducting state which is consistent with BCS-type 
mean field treatment (see supplementary). The super- 
conductivity is characterized by a finite superfluid den- 
sity Us = 1/2 which is calculated from the change of 
the ground state energy Eg upon adding a small phase 
twist 6 as [38l [40j. Moreover the finite jump for rig at 
the transition point Unnn = —0.8 (inset of Fig. [2| points 
towards a first-order transition between the valley polar- 
ized state and the spin polarized superconducting state. 
The typical momentum dependence of energy (Fig. [2jc) 
differs drastically from the 2/3 valley-polatized FQH case 
(Fig. ^) as the ground state is in the k = 0-sector with- 
out the typical quasi-degeneracy of FQH state. 

Half- filling n = PLL. We now turn to the case 
of neutral graphene (filling factor u = 0) under large 
pseudomagnetic fields generated by strain. Due to the 
electron-hole symmetry, the n = PLL is half-filled 
and there is natural a competition between valley fer- 
romagnet = Ylk k k ^P^^ ferromag- 

net ^5 = Ylk ^i? k /c tl^^ ground states, k labeling the 
Landau orbitals of the zero energy (n = 0) PLL and 
(R^L) the valleys. Similar issue of valley and spin fer- 
romagnetism in the half-filled n = Landau level has 
attracted a lot of interest for unstrained graphene under 
a real magnetic field |27[ |28l l4T|j44] . Here we revisit this 
problem in the case of a time-reversal symmetric pseu- 
domagnetic field. 

We first consider the case of pure Coulomb interaction 
{Uq — Unnn = 0). Usiug Hartrcc-Fock method ^27^128] 
we find that the valley and spin polarized states have 
the same energy when only dominant density-density 
terms are taken into account. We find that the inter- 
valley backscattering terms lift this degeneracy by favor- 
ing the valley polarization. Note that for real magnetic 
field, those backscattering terms are absent in the zero- 
energy Landau level (n = 0) due to the symmetry of 
the eigenspinors [271 EH] . Also contrary to the real mag- 
netic field Hall effect, long range Coulomb interaction 
prefers an Ising-like Z2 valley polarized state rather than 
a more general 5'/7(2)-valley-rotated state (see Supple- 
mental material). 

We now introduce on-site Hubbard interaction Uq while 
Unnn = in Eq. [2] and compute numerically the total en- 
ergy of finite size systems on a torus (Fig. [sj squares). As 
expected solely the energy of the valley polarized state is 
modified while the spin polarized state is unchanged (Fig. 
[sj horizontal dashed line) because double occupancy is 
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FIG. 3: (Color online) Neutral graphene: = 0. Upper 
panel: The Hartree-Fock energies of the Ising valley polarized 
{^v) and spin polarized ($s) states as a function of the on- 
site Hubbard coupling Uo while Unnn = (squares). We 
have also plotted the Hartree-Fock energies as a function of 
Unnn in the absence of on-site correlation, Uo = (triangles). 
Lower panel: Corresponding phase diagrams. The energy is in 
units of e^/47reao where ao is the distance between the nearest 
neighbor sites. The red dot indicates the value of the coupling 
Uo for freestanding and unstrained graphene according to Ref. 
[25] , Numerical parameters: The lattice we considered has 
96 X 96 sites. The degeneracy of the PLL orbits were Ns = 48 
while electron number is A^e = 96. 



forbidden by the Pauli principle in the fully spin polar- 
ized state. As a result, the competing valley polarized 
state (Fig. [sj empty squares) is the groundstate as long 
as the Hubbard interaction is not too repulsive {Uq < 0.5 
in units of e^/47reao), including the pure Coulomb case. 
Further increase of the on-site Hubbard interaction stabi- 
lizes the spin ferromagnet state. Using gating or different 
substrates, it could be possible to switch the groundstate 
between spin ferromagnet and valley Ising ferromagnet. 
Spin polarized STM and Kerr imaging could indeed de- 
tect these competing ground states. Here the valley Ising 
ferromagnet is an integer quantum Hall state with two 
units of quantized Hall conductance, that spontaneously 
breaks time reversal symmetry. 

In order to test the sensitivity of the phase diagram 
with respect to the details of the short range part of the 
interaction, we further consider a second model where the 
next-nearest-neighbouring coupling Unnn is varied while 
Uo = 0. Interestingly we have obtained the reverse phe- 
nomenology where repulsive Unnn tends to valley polarize 
the system (Fig. [sj triangles). 
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Conclusion. We have shown that strained graphene 
hosts various fractional topological phases which depend 
on the detailed structure of the electron-electron inter- 
actions. In current experiments on both real graphene 
[4 and artificially designed molecular graphene [5 , the 
nano-scale strained regions are small, typically of the or- 
der of the magnetic length, and they are strongly coupled 
to a metallic substrate. Future experiments on insulating 
substrates could address bigger strained regions. Never- 
theless, signatures of fractional states in restricted do- 
mains and interactions with itinerant electrons outside 
the strained region will be important topics for future 
study. The n = Landau levels are expected to be the 
best isolated, since they occur at the Dirac point, where 
the density of itinerant states is the smallest. 

The predited phases relies on the flatness of the PLL 
n = which requires spacially homogeneous strained in- 
duced magnetic fields in each valley. To this respect arti- 
ficially patterned honeycomb lattices [5, ^,'45'^ potentially 
allows for a better control upon the deformation pattern 
and therefore upon the flatness of the PLLs, in compari- 
son to the mechanical strain in real graphene. Cold atoms 
in hexagonal optical lattices [7, 8 are particularly suit- 
able to access the attractive interaction regime. Finally 
we stress that the long range part of the Coulomb in- 
teraction is always present in our calculations. This is at 
odds with current experiments [4l|5l|7l[8] but it should be 
relevant for real graphene and for future experiments on 
artificial graphenes realised in surface states lying on in- 
sulating substrates. Finally this study opens the prospect 
of discovering a series of new nontrivial topological phases 
at other fractional fillings and in higher pseudo Landau 
levels as well. 

We are grateful to J. Alicea and N. Regnault for useful 
comments. JC acknowledges support from EU/FP7 un- 
der contract TEMSSOC and from ANR through project 
2010-BLANC-041902 (ISOTOP). This work is also sup- 
ported by DOE Office of Basic Energy Sciences un- 
der grant DE-FG02-06ER46305 (DNS), the NSF grant 
DMR-0958596 for instrument (DNS), and DE-AC02- 
05CH1123 (AV) as well as the Laboratory Directed Re- 
search and Development Program of Lawrence Berkeley 
National Laboratory under US Department of Energy 
Contract No. DE-AC02-05CH11231 (PG). AV acknowl- 
edges hospitality from ICTS Bangalore where part of this 
work was completed. PG also acknowledges support from 
institute for condensed matter theory at University of 
Illinois at Urbana-Champaign. 

Note: After submission of our preprint 
|arXiv:1111.3640), a study of non zero PLLs in strained 
graphene [46] , and related work on spin fractional 
topological insulators [47 appeared. 
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SUPPLEMENTARY 
Noninteracting strained graphene 

Here we consider noninteracting spinless fermions on 
the honeycomb lattice (including of spin is straightfor- 
ward) . The triangular Bravais lattice Vmn = + ^^2 
is generated by the basis vectors: 

ai = VSaex and a.2 = -^{e^ + V^By), (3) 
and the vectors: 

^1 = ^{^^x +e^),^2 = ^{-VSej, + e^),^3 = -ae^, 

(4) 

connect any A atom to its three nearest B atoms, ao = 
0.142 nm being the length of the carbon-carbon bond. 
The area of the unit cell is Ac = 3>/3aQ/2. 

Strained induced gauge potential 

In the absence of interactions, the tight-binding Hamil- 
tonian of strained graphene (Eq. 1 in the main text) can 
be written as: 



field 6ta{rmn) of the nearest-neighbour hopping element 
between sites Vmn and Vmn + with respect to the un- 
perturbed value t. Note that the deviations Sta{rmn) are 
real quantities and must be smaller than t. 

In the absence of strain {Sta{rmn) = 0), the lattice 
Hamiltonian can be easily diagonalized and the low en- 
ergy excitations correspond to the states close to the two 
gapless Dirac points = ^(47r/3v/3<^o)6x7 where ^ = ± 
is the valley isospin. Looking for low energy effective 
theory, we expand the annihilation operators as: 

a{rmn) a+(r^n)e"'^-''""" + a_(r^n)e'^''''"", (6) 

in terms of the slowly- varying fields a^(r^n) (a similar 
equation holds for 6^(r^n) operators). Substituting into 
Eq. [5] and going from lattice to continuum fields as 

J d^y^lAc and a^(r^n) ^ \fAcH^^ ^^(rmn) 

^ ran 

(7) 

leads to the effective Hamiltonian: 

Ho = VF f d'^ Yl +4a,)^^(x), (8) 

where vp = 3at/2h is the Fermi velocity and (px^Py) = 
{—ihdx^ —ihdy) are the components of the canonical mo- 
mentum. The Pauli matrices and ay act on the lattice 
isospinors ^^(x) = (a^(x), 6^(x)). 

In the presence of a slowly varying deformation field 
^tai'Tmn) = ^^^a(x), the vallcys remain decoupled and the 
effective Hamiltonian reads: 

Ho = VF f ^'x ^l(x)(^n|a, + n^a,)^^(x), (9) 

where 11^ = P + ^ e A shows the coupling of the electronic 
charge — e with a valley- dependent gauge field A^(x) = 
^A(x) defined by: 

evF{A,{^)^iAy{^))= Y ^tai^y"^-^^' (10) 

a=l,2,3 

More specifically using Eq. [4j one obtains the gauge vec- 
tor potential in terms of the deformation field: 



^0 = ^ Y {t^^ta{^mn)){ci\rmn)K^mn^Sa)^h.C.), 

a=l,2,3 

(5) 

where second quantization operators a{rmn) and 6(r^^ + 
Sa) annihilate a fermion at A- type and B-type sites re- 
spectively. The strain is described by the deformation 



cvfA^ = --{Sti ^ dt2) ^ Sts, 
evFAy = ^{Sh-St2), (11) 
where all the position arguments x have been omitted. 
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Example of a deformation pattern and its induced 
gauge field 

We now consider a particular deformation field where 
only the bonds along ^3 = —aey are modified according 
to the pattern 



(12) 



According to Eq. 9, the corresponding strain-induced 
vector potential is the familiar vector potential of the 
Landau gauge A = —Bye^^ describing here a uniform 
magnetic field B = Bbz in the valley ^ = +, and the 
opposite field in valley ^ = — . In this gauge, the natural 
geometry is a rectangular one with linear sizes and Ly. 
Due to translational invariance, the system can be infinite 
in the x-direction but the smooth deformation condition 
{Sts <C t) brings a limitation on the transverse size Ly 
because Sts is growing linearly long the ^/-direction. As- 
suming 6ts{y = 0) = 0, then dts{Ly) = evpBLy cannot 
exceed a reasonable fraction oit = 2hvF/^CL which leads 
to the condition: 



^1 

\B\a 



(13) 



Typically for a magnetic length Ib = \/h/eB 10 nm 
and lattice constant ao — 0.1 nm, the ribbon width can- 
not exceed 500 nm. The strength of the effective mag- 
netic field B is proportional to the gradient of the hop- 
ping amplitude deformation. For a similar global defor- 
mation over the whole sample, a narrow ribbon hosts a 
stronger magnetic field than a broader ribbon. Note that 
many other deformation fields lead to the same gauge 
vector potential, including {Sti, St2, Sts) = evFBy{2^ 0, 0) 
or {6ti, St2, Sts) = evpByiO, 2,0). 

Single electron wavefunctions 

We derive here the wavefunctions for noninteracting 
Dirac fermions under a strong pseudo-magnetic field 
(^^e^), or more precisely the valley-dependent gauge po- 
tential A^(x) = ^A(x), which are both opposite fields in 
the valleys ^ = ±. We choose B positive for definite- 
ness, and denotes Ib = \/h/eB the magnetic length. For 
each valley, we consider the first quantized Hamiltonian 
corresponding to Eq. [9j namely: 



(14) 



where the components of the gauge-independent momen- 
tum = p -|- ^eA do not commute due to the presence 
of the pseudo-magnetic field. Unlike the real magnetic 
field case, the sign of the commutator: 



[n|,n«] 



(15) 



depends on the valley index ^. Hence the ladder opera- 
tors are defined as 

Ib 



V2h 



(n« -icn«), 



(16) 



in order to enforce the proper commutation relation 



^tl _ rj^j^g Hamiltonian can be written as: 







(17) 



We now focus on the zero energy Landau level. The cor- 
responding wave function in the ^-valley is (0,^^) with 
a^v^ = 0. Hence in the zero energy level, single electron 
wavefunctions are finite only on one triangular sublat- 
tice, here the B-atoms sublattice (since we have chosen 
the field strength B to be positive). This is a general 
property valid for any strain-induced gauge field on the 
graphene lattice. 

Now we give explicitly the wavefunctions for the Lan- 
dau gauge field A = —Bye^^. Then the equation a^v^ = 
reads: 



v^{x,y) = 0, 



and substituting v^{x,y) = fdv)^ we get 



dy 



^B 



kiy) = 0- 



Finally the normalized wavefunction is: 



v^{x,y) 



1 



V^'^'LJb 



(18) 



(19) 



(20) 



Those wavefunctions correspond to the continuum 
model. For the numerical calculations performed on the 
honeycomb lattice and discussed in the main text, the 
wavefunctions are different from Eq. |20] and have to be 
determined numerically using lattice model with pseudo 
magnetic field. Nevertheless for the n = 0-Landau level, 
there is a simple procedure to get the wavefunctions in 
the pseudo gauge field ^A from the ones in the usual 
valley-independent gauge field A. As said before for 
the real magnetic field, the (^ = +)-valley wavefunctions 
have predominant weight on B atoms, and (^ = — )-valley 
wavefunctions on A atoms. To get the corresponding 
wavefunctions in a pseudo magnetic field, one has simply 
to i) keep the former (^ = +) -valley wavefunctions on B 
atoms without any change, ii) remove the wavefunctions 
on A sites and finally iii) add complex conjugates of B 
sites wavefunctions. 



Effect of interactions 

Here we will consider the projection of density-density 
interactions of the form Hj = ^, V{r — r')n(r)n(r') 
into the n = PLL which then would have the form: 
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Pl,P2,P3,P4,cr,cr' 

(21) 

where Cp^ is the creation operator of fermion in state 
P = k} which is in vahey ^ and Landau orbital k (see 
eqn. |20|. 

Using the continuum model derived in the last section 
we get the the projected potential: 



Vp„PM = (^) / drdr'Vir - r')(-i)«3+«-«-«^ 
g-ifj[(fei«'-y)'+(fe2('-y')'+(fc3('-!/')"+(M"-y)"] 

(22) 

where g = (gf^) ^ • We can now consider different inter- 
action potentials. 

If the interaction is smooth the dominant interaction 
term has = ^4 and ^2 = ^3- 

There are two polarized states which will be particu- 
larly of our interest, valley {ipv) = Il/ecJ ^ ^ ^ or spin 

\ips) = Ilkc\ ^c^_i ^ polarized states. The energy of 
these states have the general form of: 



= Yl ^{P1,1},{P2,1},{P2,1},{P1,1} 

{pi,l},{p2,l},cr,cr' 

- '^^{P1,1},{P2,1},{P1,1},{P2,1}^^,^^ 

(23) 

and 

^S= Yl ^{Pl,a},{P2,42},{p2,C2},{pi,6} 

{Pl,'^l},{P2,^2} 

~ '^{pi,'ei},{p2,<e2},{pi,'ei},{P2,^2}%,^2- 

(24) 



Coulomb interaction 

Putting the Coulomb interaction V{Ti — 

rj)n{ri)n{rj) = e^n(ri)n(r^)/47re|ri - rj| in projec- 
tion 122 



Pl,P2,P3,P4 



5V 



drdr' 



r — r 



For neutral graphene the n = PLL is at half filling. 
The interactions naturally prefers the polarized state of 



valley or spin degrees of freedom. Using |24] the energy of 
spin polarized state is given by: 



2 2 \ ^ 

Es=Y. (^) I d{x-x') dydy' 



Lx I 



2 (2e-^ - cos [(pi - p2){x - x')] 



(26) 



Interestingly the valley polarized state has the same en- 
ergy: Ey = If we include the fast oscillation which 
we ignored before (e*^ ['^^i"^4)ic+(^2-^3)a^'] ^.j^^ valley po- 
larized state does not change since ^1 = ^2 = ^3 = ^4^ 
but the energy of spin polarized state increases, so the 
valley polarized state will be the ground state of neutral 
graphene if the long range Coulomb is the dominating 
interaction. 

Notice that we considered a valley polarized state in a 
single valley. If we consider an arbitrary rotated valley 
polarized state in a state of the form: 



(27) 

the Hartree-Fock energy have the form: 



(28) 



P1/P2 



g - i [<^i ( /e 1 - /e4 ) +<^2 ( ^2 - /C3 ) a^' ] 



(25) 



where Tp^^p^ is positive for the Coulomb potential. The 
term added is then always positive and the minimum of 
energy is for ii+ = or = for all p. The the preferred 
state would be an Ising valley polarized state. 



Short range Hubbard interactions and 
superconductivity 



Although all types of density-density interactions could 
be treated using our projection scheme, here we only 
present the details for the on-site and next nearest neigh- 
bor interactions which we also studied numerically and 
observed the somehow unexpected spin polarized super- 
conductivity. 

The on-site interaction has the form j drdr'VdiY — y') 
where as next-nearest neighbor interaction has the form 
/ dvdv'V YH^^ (5(r + - v') where = (X^, Yi) are the 
vectors connecting each site to its six next nearest neigh- 
bors. This next nearest neighbour interaction projected 
into the n = PLL has the form: 
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The corresponding mean-field equations for the super- 
conducting gap for the spin singlet state reads: 



dy 



^(6-6). 

(29) 



Here we have 6 - 6 = ^4 - 6 and ^ipi + ^2P2 - ^sPs - 
^4P4 = 0. Notice that the interaction is not smooth in 
the lattice scale so we should keep the oscillatory term 
m2 - 6) = Ell e^tf fe-«3)x. ^ 2 COS [f (6 - 6)] + 
4 cos [f (6-6)]. 

For on-site interaction we have similar form with 

f^(C2-6) = i. 

With this form of the projected interaction potential 
we can readily compare the energy of spin and valley 
polarized states. Using the expressions in 23 and 24 we 
get: 



Ev =V 



Lx PJ"^' ^ 

Pl,P2 



(30) 



and 



Pl,P2 



(31) 



For the on-site interaction with ^{^2—^1) = 1 the 
spin polarized state has no energy-gain where as valley 
polarized state energy increases with on site repulsion. 

For the next nearest neighbour interaction ^^(^2—^1) = 
(^(0)(|^^^,^2 ~ 1) < ^^i,'^2- valley polarized state will 
be stabilized with next nearest neighbor interactions. 

With the projected interaction of the general form 

^Pl,P2,P3,P4,a,a' ^Pl,P2,P3,P4 4,Pi4',P2^^''^3Ca,P4 ^C CaU 

also compare the mean-field energy of different supercon- 
ducting states. We compare the energy of two supercon- 
ducting states which are valley triplet, spin singlet: 



l^^Singlet) =Up^=^^^p^y^^_^p^y{Ui^ViC^p^^^C^p^^^)\0), (32) 

and spin triplet, valley singlet: 

{'^Triplet) = Up^ = {^^p^y^^=^^^{Ui + ^icj.^ a) 1^) ' (^^) 

where P and a are the time reversal of P and a respec- 
tively. 



Ap = — Vppp/p/Up'Vp' 
and for the spin triplet state reads: 



Ap = — Vp p pi piUp'Vp' 



(34) 



(35) 



where 



dy e 



2 cos 



47r 



(6-6) 



4 cos 



27T 



(6 - 6) 



In the spin-singlet state Vj 



p,p,P',P' 



changes sign be- 
tween the states with ^ = and 6" =6 where as for 
the spin triplet for all the term in the mean field equation 
have 6=6 So we get larger gap and so smaller mean- 
field energy for the spin triplet superconducting state. 



The robustness of the fractional valley Hall insulator 
and its evolution with tunning Unnn 

The fractional valley Hall insulator is characterized by 
a large ground state degeneracy. In a finite-size system, 
due to the coupling between different states, one usually 
sees a quasi-degeneracy with a finite splitting between 
these states in the ground state manifold. As shown the 
main text, around U^^^ = — 0.6 without the interaction 
between the electrons in the same valley U^^^ = 0, the 
nine states at the right quantum number sectors are in- 
deed have much lower energy than other excited states. 
However, the splitting between these states are close to 
the finite gap between these states and other excited 
states. While the quantized nonzero total Chern number 
of the ground state manifold indicates the obtained state 
is indeed a fractional valley Hall insulator, a clear-cut ev- 
idence of nine fold topological degeneracy is still absent 
and it is difficult to predict the fate of the state as system 
becomes very large. Here we address this issue through 
tuning the system deep into the topological phase. In- 
deed this can be achieved by increasing the correlations 
between the electrons in the same valley {U^^^). 

As shown in Fig. [4j with the turn on of positive U^nn^ 
the energy gap between the ground state manifold and 
other excited states becomes very robust and much larger 
than the splitting of the energy of the ground state man- 
ifold. There is no phase transition as Unnn continuously 
increases, so the observed state is indeed the same phase 
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FTI, U°Pnnn=-0.6, Nl,T=Nr,T=4 




-1 -0.5 O 0.5 1 1.5 2 2.5 3 

Hubbard Unnn^ 



as the fractional valley Hall insulator at the decoupled 
limit (strong U^^^ limit). In that limit, the nine fold de- 
generacy is exact and spinless electrons in different val- 
leys are contributing ±1/3 quantized Hall conductances. 
We further perform the flux inserting measurement, we 
show in Fig. for U^^^ = 1, the three lowest energy 
states in the momentum sector k = evolve into other 
states in the ground state manifold and they evolve back 
to themselves after three periods of boundary phase in- 
sertion. The energy gap between these states and other 
excited states remain robust as illustrated in Fig. ^p. 
We further perform valley-dependent Chern number cal- 
culation, and find a total Chern number quantized to 6 
for all nine levels, characterizing the 2/3 fractionalized 
valley spin-Hall effect. Remarkably, this phase persists 
in a wide range of U^^^ > —0.2 including the simple case 
where this interaction is turned off = as shown 

in the main text part of the paper. 




FIG. 4: (Color online) (a) Fractional topological insulating 
phase for spinless electrons at filling ly — —2 + 2/3. The 
LL degeneracy is Ns — 12 while there are totally Ne = 8 
(Nl — Nr — 4) electrons with polarized spin on a lattice with 
24 * 24 sites. We demonstrate an emergent symmetry where 
energies of nine states from k — 0, 27r/3 and 47r/3 sectors 
become near degeneracy at large Unnn limit. The onset of the 
fractionalized phase with nine- fold near degenerating GSM is 
identified at Unnn ^ —0.2. (b) The evolution and the robust 
of the gap during the change of the boundary phase. 



